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[i]  Spatial  and  temporal  variability  of  the  impact  of  air-sea  stratification  on  the 
differences  between  satellite-derived  10  m  equivalent  neutral  wind  speeds  and  stability- 
dependent  (e.g.,  in  situ)  10  m  wind  speeds  are  quantitatively  examined  over  the  global 
ocean.  The  influences  of  stability  are  compared  with  three  air-sea  flux  algorithms,  Coupled 
Ocean- Atmosphere  Response  Experiment  (version  3.0),  Bourassa- Vincent- Wood,  and 
Liu-Katsaros-Businger.  Analyses  are  first  presented  at  many  individual  buoy  locations  and 
then  are  extended  to  the  global  ocean  with  the  use  of  rain-free  wind  measurements 
from  the  SeaWinds  scatterometer  on  the  QuikSCAT  satellite,  gridded  at  a  resolution  of 
0.25°  x  0.25°.  Overall,  stability-dependent  winds  are  found  to  be  weaker  than  equivalent 
neutral  winds  by  0.2  ms-1  on  the  basis  of  7619  monthly  mean  values  from 
208  buoys  during  2000-2005.  Differences  based  on  hourly  winds  can  be  as  large  as 
±0.5  ms-1.  Results  remain  robust  regardless  of  which  air-sea  flux  algorithm  is  used. 

Monthly  rain-free  gridded  QuikSCAT  measurements,  combined  with  atmospheric  stability 
determined  using  near-surface  variables  from  the  European  Centre  for  Medium-Range 
Weather  Forecasts  40-year  reanalysis,  demonstrate  the  effects  of  stratification  on  the  10  m 
winds  globally.  Differences  in  stability-dependent  and  neutral  winds  are  substantially 
nonsymmetrical  and  reveal  locations  where  the  former  is  stronger  than  the  latter.  These 
differences  may  cause  physically  significant  biases  in  air-sea  fluxes  if  they  are  not 
properly  considered,  especially  near  the  Kuroshio  and  Gulf  Stream  current  systems. 

Citation:  Kara,  A.  B.,  A.  J.  Wallcraft,  and  M.  A.  Bourassa  (2008),  Air-sea  stability  effects  on  the  10  m  winds  over  the  global  ocean: 
Evaluations  of  air-sea  flux  algorithms,  J.  Geophys.  Res.,  113 ,  C04009,  doi:  10.1 029/2007 JC004324. 


1.  Introduction 

[2]  Accurate  wind  speeds  are  essential  for  reliable  com¬ 
putations  of  surface  heat  and  momentum  fluxes  (e.g.,  wind 
stress  and  sensible  and  latent  heat  fluxes)  over  the  global 
ocean.  The  sampling  density  (particularly  for  data-sparse 
regions)  and  accuracy  make  satellite  winds  desirable  data 
for  many  related  applications  such  as  coastal  upwelling, 
oceanic/atmospheric  coupling  associated  with  both  tropical 
instability  wave  and  ocean  fronts  [Chelton  et  a/.,  2004], 
ocean  currents  [Lagerloef  et  al. ,  2003],  and  detection  of 
tropical  disturbances  [Gierach  et  al. ,  2007],  to  list  a  small 
sample  of  applications. 

[3]  Quantifying  the  effects  of  atmospheric  stability  on 
satellite-derived  wind  speed  (called  equivalent  neutral  wind 
speed)  is  a  critical  issue  for  various  climate  applications, 
where  biases  related  to  such  processes  would  be  of  interest 
in  flux  calculations,  product  intercomparisons,  and  ocean 
model  simulations.  Of  particular  interest  are  scatterometers, 
spacebome  radars  that  infer  surface  winds  from  the  rough- 
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ness  of  the  ocean  surface  [Liu,  2002].  They  respond  to 
changes  in  the  water  surface  rather  than  responding  directly 
to  wind  speed. 

[4]  Surface  roughness  is  more  closely  related  to  wind 
stress  than  to  wind  speed  [Bourassa,  2006].  While  this 
would  seem  ideal  for  ocean  forcing,  the  scatterometer  time 
series  has  existed  for  a  relatively  short  time,  and  most  ocean 
models  are  forced  by  stress  determined  from  wind  speed 
[e.g.,  Barron  et  al.,  2006].  Satellite  winds  (e.g.,  from 
scatterometers  and  radiometers)  are  well  suited  for  such 
applications,  provided  that  considerations  such  as  ocean 
currents  [Kara  et  al.,  2007a]  and  atmospheric  stability 
[Kara  et  al.,  2005]  are  properly  considered  in  the  conver¬ 
sion  of  satellite  winds  to  stresses.  Satellite  winds  are 
calibrated  such  that  these  considerations  have  already  been 
accounted  for:  further  adjustments  result  in  regional  biases 
in  surface  forcing.  Specifically,  scatterometers  and  radio¬ 
meters  are  calibrated  to  equivalent  neutral  wind  speeds  at 
10  m  above  the  ocean  surface  [e.g.,  Meissner  et  al.,  2001] 
rather  than  to  wind  speeds  including  effects  of  stratification. 
Hereinafter,  the  latter  will  be  referred  to  as  stability-depen¬ 
dent  winds. 

[5]  The  conversion  of  equivalent  neutral  wind  speeds  to 
stability-dependent  ones  is  complicated  by  several  defini¬ 
tions  of  the  neutral  wind  as  well  as  conflicting  assumptions 
used  to  determine  an  equivalent  neutral  wind.  In  general, 
there  are  two  definitions.  In  the  most  common  definition. 
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Table  1.  List  of  Algorithms  Used  in  This  Paper® 


Abbreviation 

Air-Sca  Flux  Algorithm 

Reference 

COARE 

Coupled  Occan-Atmosphcrc 
Response  Experiment 

Fairall  et  al.  [2003] 

LKB 

Liu-Katsaros-Businger  model 

Liu  and  Tang  [1996] 

BVW 

Bourassa-Vincent-Wood  model 

Bourassa  et  al.  [  1 999] 

BVWN 

Bourassa-Vincent-Wood 
neutral  model 

Bourassa  et  al.  [1999] 

LOG 

Logarithmically  varying  profile 

Peixoto  and  Oort  [  1 992] 

“Abbreviations  used  throughout  the  paper  are  given.  A  brief  description 
for  each  algorithm  is  provided  in  the  text.  For  detailed  descriptions,  see 
references.  The  modem  version  of  the  LKB  model  [Liu  and  Tang ,  1996] 
differs  from  the  original  algorithm  in  that  it  increases  the  momentum 
roughness  length  by  adding  roughness  due  to  gravity  wave  via  a  Chamock 
parameterization. 


equivalent  neutral  wind  speed  is  the  mean  wind  speed  that 
would  be  observed  if  there  was  neutral  atmospheric  strati¬ 
fication  [Geemaert  and  Katsaros,  1986].  This  definition  is 
useful  for  determining  neutral  drag  coefficients  and  rough¬ 
ness  lengths.  However,  this  definition  requires  assumptions 
on  what  (if  anything)  is  held  constant,  which  may  not  be 
appropriate  globally.  The  second  definition  [Ross  et  al , 
1985;  Liu  and  Tang ,  1996]  of  equivalent  neutral  wind  is  the 
wind  speed  calculated  by  using  the  stress  and  roughness 
length  consistent  with  the  observed  atmospheric  stratifica¬ 
tion  but  setting  the  atmospheric  stratification  term  in  the 
modified  log-wind  profile  equal  to  zero.  This  definition  may 
be  more  consistent  with  scatterometry,  which  measures 
surface  stress  [Weissman  et  al .,  1994].  Unfortunately,  accu¬ 
rate  measurements  of  near-surface  stress  over  water  are 
extremely  sparse  relative  to  wind  speed.  Another  issue  to 
consider  is  that  scatterometer  equivalent  neutral  winds  are 
relative  to  current  and  are  further  modified  by  surface  wave 
motion  [ Bourassa ,  2006]  from  scatterometer  winds.  The 
considerations  of  wave  motion  and  currents  are  physically 
important  in  the  calculation  of  surface  turbulent  fluxes  [e.g., 
Kara  et  al .,  2007a].  Herein,  we  focus  on  the  impact  of 
atmospheric  stability. 

[6]  There  are  earlier  studies  related  to  satellite  calibration/ 
validation,  which  adjusted  in  situ  observations  to  be  phys¬ 
ically  consistent  with  equivalent  neutral  winds.  For  exam¬ 
ple,  Mears  et  al.  [2001]  used  Special  Sensor  Microwave 
Imager  winds  during  1987-1997,  and  Ebuchi  et  al.  [2002] 
analyzed  QuikSCAT  winds  during  1999-2000  at  many 
individual  buoy  locations.  Bourassa  et  al.  [2003]  evaluated 
SeaWinds  data  with  research  vessels.  Each  of  these  earlier 
studies  adjusted  the  in  situ  observations  to  equivalent 
neutral  winds.  Mears  et  al.  [2001]  found  that  the  distribu¬ 
tion  of  adjustments  had  little  dependence  on  the  flux 
algorithm  used  for  adjustment.  However,  they  expressed 
this  distribution  through  means  and  standard  deviations.  We 
will  show  that  the  distribution  of  adjustments  is  substan¬ 
tially  non-Gaussian  and  that  there  can  be  large  changes  in 
the  distribution  if  an  overly  simple  adjustment  algorithm  is 
applied. 

[7]  In  this  paper,  we  focus  on  the  distribution  of  such 
differences  and  examine  the  regional  variability  on  a  fine 
(relative  to  most  global  ocean  models)  grid  resolution 
(0.25°  x  0.25°)  over  the  global  ocean  during  2000-2005. 
This  is  accomplished  by  using  a  variety  of  air-sea  flux 
algorithms.  We  quantify  the  difference  between  equivalent 


neutral  winds  and  those  including  influences  of  stratification 
to  provide  a  comprehensive  analysis  over  the  global  ocean. 
This  is  the  first  study  demonstrating  how  the  differences 
vary  regionally  and  globally. 

2.  Air-Sea  Flux  Algorithms 

[8]  The  World  Meteorological  Organization  has  defined 
10  m  as  the  reference  height  for  near-surface  wind  obser¬ 
vations.  We  will  use  various  algorithms  to  adjust  wind 
speed  at  an  observed  height  (e.g.,  4  m  above  the  surface) 
to  the  standard  height  of  10  m.  Each  algorithm  will  be 
applied  to  observed  winds,  some  of  which  include  air-sea 
stratification  effects  and  some  of  which  do  not. 

[9]  A  list  of  algorithms  used  for  adjusting  observed  winds 
to  1 0  m  above  the  ocean  surface  and  calculating  equivalent 
neutral  winds  and  stability-dependent  winds  is  provided  in 
Table  1  along  with  their  abbreviations.  The  reason  for  using 
different  algorithms  is  to  investigate  whether  or  not  the 
resulting  adjusted  winds  deviate  significantly.  These  algo¬ 
rithms  have  been  well  described  in  the  literature.  Thus,  only 
a  brief  description  for  each  is  provided.  Further  details  can 
be  found  in  references  (see  Table  1). 

[10]  The  Coupled  Ocean- Atmosphere  Response  Experi¬ 
ment  (CO ARE  version  3.0)  and  Bourassa- Vincent- Wood 
(BVW)  algorithms  will  be  used  for  obtaining  height-adjusted 
stability-dependent  10  m  winds  (i.e.,  including  stratification 
effects  at  the  air-sea  interface).  The  Liu-Katsaros-Businger 
(LKB),  Bourassa-Vincent-Wood  neutral  (BVWN),  and  log¬ 
arithmically  varying  profile  (LOG)  algorithms  will  be  used 
for  obtaining  the  10  m  equivalent  neutral  winds.  Unlike  the 
COARE  algorithm,  BVW  has  an  option  of  calculating  both 
equivalent  neutral  winds  and  winds  based  on  the  air-sea 
stability.  All  algorithms  have  advantages  and  disadvantages 
over  each  other  because  of  their  unique  parameterizations 
[e.g.,  Brunke  et  al. ,  2003],  whose  details  are  not  given  here. 

[11]  The  COARE  algorithm  attempts  to  take  full  account 
of  atmospheric  stability  in  the  near-surface  layer.  This  is 
achieved  by  iteratively  solving  for  the  stability  of  the 
atmosphere.  The  algorithm  uses  the  stable  profile  functions 
from  Beljaars  and  Holtslag  [1991]  and  the  convective 
portion  of  the  scalar  profile  function  described  by  Grachev 
et  al.  [2000].  The  threshold  value  z/L  =  1,  where  z  is  the 
height  above  the  ocean  surface  and  L  is  the  Monin-Obukhov 
length,  for  very  stable  cases,  is  not  a  strict  limit  in  the 
algorithm.  The  use  of  the  critical  bulk  Richardson  number 
is  more  suitable  to  setting  a  very  stable  limit.  A  bulk 
Richardson  number  of  about  0.2  is  taken  as  the  critical  value 
[Grachev  et  al .,  2005].  There  is  no  explicit  option  of  the 
equivalent  neutral  winds  in  the  algorithm. 

[12]  The  BVW  model  is  a  fully  coupled  flux  and  ocean 
state  model.  Coupling  allows  the  ocean  state  to  respond  to 
atmospheric  stability  and  vice  versa.  The  profile  stability 
functions  in  the  algorithm  were  obtained  from  Beljaars  and 
Holtslag  [  1 99 1  ]  for  stable  conditions  and  from  Benoit  [  1 977] 
for  unstable  conditions.  Typically,  the  parameterizations  in 
BVW  relate  surface  roughness  lengths,  and  hence  the 
exchange  coefficients,  to  various  aspects  of  sea  state,  swell, 
gravity  waves,  and  capillary  waves.  It  can  provide  both 
stability-dependent  and  equivalent  neutral  winds. 

[13]  LKB  is  also  based  on  the  similarity  functions  (i.e., 
nondimensional  flux-profile  relations)  in  the  atmospheric 
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Figure  1.  (a)  Typical  variations  of  (top)  actual,  i.e., 
stability-dependent,  and  (middle)  equivalent  neutral  winds 
with  height  along  with  (bottom)  the  differences  between  the 
two.  The  black  lines  for  stability-dependent  and  neutral 
winds  are  similar  because  they  represent  conditions  near 
neutral  atmospheric  stratification,  (b)  Comparison  of 
stability-dependent  and  neutral  wind  speeds  for  (top) 
unstable,  (middle)  neutral,  and  (bottom)  stable  conditions. 
Stable  (unstable)  stratification  corresponds  to  positive 
(negative)  air-sea  temperature,  and  it  reduces  (increases) 
mixing  through  surface  stress. 


surface  layer.  As  in  COARE  and  BVW,  bulk  parameter- 
izations  are  used  for  determining  flux-profile  relationships 
of  heat  and  moisture  within  the  surface  layer.  The  interfacial 
effects  within  the  molecular  sublayer  at  the  sea  surface  are 
included  because  it  is  assumed  to  be  the  most  limiting 
(slowest)  process.  The  wind  profile  in  the  sublayer  profile  is 
given  an  exponential  shape.  Theoretically,  the  validity  of 
LKB  depends  on  the  validity  of  the  modeled  surface 
renewal  and  similarity  theory.  For  example,  when  turbu¬ 
lence  is  suppressed  by  stable  density  stratification  (bulk 
Richardson  number  exceeds  a  critical  value),  the  results 
may  not  be  valid  and  the  algorithm  may  fail  to  converge. 

[14]  The  LOG  approach  is  the  simplest  one  that  is  used 
for  adjusting  winds.  Similar  to  LKB,  it  works  for  the 
equivalent  neutral  winds.  Assuming  a  logarithmically  vary¬ 


ing  wind  profile,  the  corrected  wind  speed  at  a  height  z  is 
given  by  ULOG(z)  =  \n(z/z0)/[\n(zjz0)  (J(Zm)].  Here,  za  is 
the  roughness  length.  The  typical  oceanic  value  for  z0 
is  1.52  x  10~4  m  [Peixoto  and  Oort ,  1992]. 

3.  Definition  of  the  Equivalent  Neutral  Wind 

[15]  The  wind  speed  dependence  on  height  is  typically 
expressed  as  follows: 

U(z)  -Us  =  (M*/*)[ln(z/z0)]  +  (m*/*)<I>(z,z0,I),  (1) 

where  term  1  (U(z)  -  Us)  denotes  the  wind  speed  relative  to 
the  surface  speed  Us  (e.g.,  the  surface  current)  as  a  function 
of  height  z.  In  term  2  ((u*/k)[\n(z/z0)])>  the  friction  velocity 
(u*)  is  the  square  root  of  the  kinematic  stress,  and  the  von 
Karman  constant  (k)  has  a  value  of  0.4.  Term  2  is  the 
adjustment  for  height  that  does  not  directly  consider 
atmospheric  stability;  however,  atmospheric  stability  enters 
indirectly  through  the  friction  velocity  and  roughness  length. 
The  quantity  4>(z,  z0)  in  term  3  ((uJk)${zy  z0,  L))  denotes  the 
direct  modification  due  to  atmospheric  stratification  based 
on  the  Obukhov  scale  length  ( L ).  Both  u*  andzG  are  functions 
of  Us ,  atmospheric  stratification,  and  ocean  state.  When  the 
atmospheric  stratification  is  neutral  (i.e.,  z/L  =  0),  there  is  no 
stratification,  and  therefore,  the  stability  term  is  set  to  zero. 
Equivalent  neutral  winds  are  calculated  without  term  3; 
however,  they  use  roughness  length  and  friction  velocity  that 
are  first  determined  with  term  3  included. 

[16]  Following  (1),  we  illustrate  the  influence  of  air-sea 
temperature  contrasts  on  the  wind  profile  for  both  winds, 
including  the  effects  of  stability  and  the  corresponding 
equivalent  neutral  winds  (Figure  la).  In  the  example,  the 
wind  speed  is  set  to  6  m  s_l  at  a  height  of  6  m.  The  air-sea 
mixing  ratio  difference  is  taken  as  3  g  kg-1  over  the  same 
height  range,  causing  the  profile  with  a  zero  temperature 
difference  to  be  slightly  unstable. 

[17]  Equivalent  neutral  wind  profiles  given  in  Figure  la 
are  determined  with  the  definition  of  Liu  and  Tang  [1996]. 
Clearly,  there  are  considerable  differences  between  the  two 
sets  of  profiles.  The  wind  profile  is  equal  to  the  equivalent 
neutral  wind  profile  only  for  the  case  of  neutral  stratifica¬ 
tion.  Typically,  these  differences  are  <0.5  m  s_I;  however, 
they  can  be  very  large  for  atmospheric  stratifications  that  are 
far  from  neutral,  as  evident  from  Figure  la  (bottom). 

[is]  For  stable  atmospheric  stratification  (green  and  blue 
lines),  equivalent  neutral  winds  are  relatively  weaker  than 
stability-dependent  winds,  and  for  unstable  atmospheric 
stratification  (red  and  cyan  lines),  equivalent  neutral  winds 
are  stronger  than  stability-dependent  winds.  Near  the  sur¬ 
face  (i.e.,  below  the  top  of  the  boundary  layer  with  increas¬ 
ing  differences  closer  to  the  surface),  the  wind  speed  is 
typically  weaker  when  the  air  is  stable  and  stronger  when 
the  air  is  unstable  (Figure  lb). 

4.  Buoy  Data  and  Wind  Speed  Adjustment 
to  10  m 

[19]  We  use  buoy  data  from  three  sources:  (1)  the 
Tropical  Atmosphere-Ocean  (TAO)ATriangle  Trans-Ocean 
Buoy  Network  (TRITON)  array  [McPhaden  et  al. ,  1998], 
(2)  the  Pilot  Research  Moored  Array  in  the  Tropical  Atlantic 
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Figure  2.  Locations  of  TAO  (crosses),  PI  RATA  (pluses), 
and  NDBC  (circles)  buoys  used  in  the  analyses  in  this  paper. 
The  regions  where  ice  exists  are  shown  in  gray.  The  wide 
distribution  of  the  buoys  and  their  records  provides  an 
ample  source  of  data  from  which  to  test  algorithms  used  in 
this  paper  in  both  open  ocean  (TAO  and  PIRATA  buoys  and 
some  NDBC  buoys)  and  coastal  regions  (most  NDBC 
buoys)  during  long  time  periods. 

(PIRATA)  [Servain  et  al.y  1998],  and  (3)  the  National  Data 
Buoy  Center  (NDBC)  database,  which  is  available  from  the 
National  Oceanographic  Data  Center  (http://www.nodc. 
noaa.gov/BUO  Y  /buoy.html). 

[20]  Locations  of  buoys  are  shown  in  Figure  2.  The  heights 
of  sensors  at  NDBC  locations  are  generally  the  same,  with 
values  of  4  m  for  air  temperature  and  humidity  and  5  m  for 
wind  speed  measurements.  However,  there  are  a  few  excep¬ 
tions.  For  example,  some  NDBC  buoys  measure  air  temper¬ 
ature  and  wind  speed  at  3.2,  10.0,  12.3,  and  13.8  m  above  the 
ocean  surface.  Both  TAO  and  PIRATA  buoys  have  identical 
measurement  heights:  air  temperature  and  relative  humidity 
are  measured  at  3  m  above  the  ocean  surface,  and  wind  speed 
is  measured  at  4  m  above  the  ocean  surface.  TRITON  buoys 
in  the  Pacific  west  of  160°E  measure  winds  at  3.5  m  and 
measure  relative  humidity  and  air  temperature  at  2.2  m. 

[21]  Measurements  of  ocean  surface  temperature,  air 
temperature,  and  relative  humidity  (or  specific  humidity) 
are  used  for  adjusting  winds  to  10  m  from  their  original 
heights  and  for  calculating  air-sea  stability.  Each  buoy 
provides  a  time  series  of  near-surface  atmospheric  variables, 
with  the  time  period  sampled  in  each  record  varying 
according  to  buoy  deployment  duration  and  sensor  opera¬ 
tion.  We  obtain  time  series  of  these  historical  atmospheric 
variables  at  hourly  resolution  from  all  available  buoys  from 
2000  through  2005.  There  are  78  TAO,  117  NDBC,  and 
13  PIRATA  buoys  used  in  this  study. 

[22]  For  consistency  at  all  TAO,  NDBC,  and  PIRATA 
buoys,  mixing  ratio  values  for  air  (using  air  temperature) 
and  ocean  (using  ocean  surface  temperature)  are  calculated 
using  saturated  vapor  pressure.  Air  and  ocean  surface 
temperatures  and  ocean  level  pressures  available  at  each 
location  are  used  for  computations.  To  calculate  mixing 
ratio,  dew  point  temperature  is  available  from  NDBC  buoys, 
while  TAO  and  PIRATA  buoys  provide  relative  humidity. 
TAO  and  PIRATA  buoys  do  not  generally  measure  atmo¬ 
spheric  pressure;  thus  a  constant  value  of  1 0 1 3  mbar  is  used 
for  density  calculations. 

[23]  Positions  of  ocean  buoys  can  change  by  up  to  «3  km 
over  the  course  of  a  few  days  to  a  week  depending  on  the 


current  regime.  Since  each  mooring  moves  in  time  and 
space  from  its  deployment  position,  we  calculated  average 
position  on  the  basis  of  the  historical  latitude  and  longitude 
data  for  each  buoy.  For  ease  of  notation,  nearest  integer 
values  of  average  latitude  and  longitude  are  used  for  each 
buoy  throughout  the  text. 

[24]  We  first  investigate  how  much  wind  speed  changes 
when  adjusting  it  to  10  m  from  a  particular  measurement 
height.  All  algorithms  listed  in  Table  1  are  modified  (if 
needed)  such  that  a  measurement  height  for  each  buoy  can  be 
input  for  adjusting  winds  to  10  m  as  examined  here  or 
calculating  equivalent  neutral  winds.  Computations  are  per¬ 
formed  on  the  basis  of  the  available  hourly  buoy  measure¬ 
ments,  involving  the  above-mentioned  near-surface 
atmospheric  variables.  The  resulting  winds  are  then  averaged 
over  a  day  and  a  month.  We  do  not  use  buoy  measurements 
for  periods  where  any  meteorological  parameter  required  by 
the  algorithms  is  missing,  with  the  exception  of  pressure, 
which  is  set  to  1013  mbar.  Additionally,  no  averaging  or 
smoothing  is  applied  to  the  original  parameters  prior  to  using 
them  in  the  algorithms.  Daily  and  monthly  means  are  then 
constructed  to  reduce  data  volume  at  all  buoy  locations. 

[25]  An  example  of  the  adjustment  procedure  is  demon¬ 
strated  at  a  TAO  buoy,  located  at  the  equator  (0°N,  155°W), 
in  2000  (Figure  3a).  Wind  speed  measurements  are  made  at 
a  height  of  4  m  at  this  particular  buoy  location.  As 
mentioned  earlier,  near-surface  atmospheric  variables  are 
input  to  the  COARE  and  BVW  algorithms  to  adjust  winds 
to  a  height  of  10  m.  The  change  in  wind  speed  due  to  height 
adjustment  is  proportional  to  the  observed  wind  speed. 
Therefore,  the  impact  of  adjustment  for  weaker  winds 
(e.g.,  3  m  s-1)  is  typically  small  in  comparison  to  that  for 
stronger  winds  (e.g.,  8  m  s“!).  For  example,  in  the  former 
(latter)  case  the  10  m  wind  is  «0.3  m  s_I  (w0.8  m  s_1) 
stronger  than  the  4  m  wind.  Wind  speed  at  the  measurement 
height  of  4  m  has  an  annual  mean  value  of  6.3  ms-1,  and 
wind  speed  adjusted  to  10  m  using  either  COARE  or  BVW 
has  an  identical  mean  value  of  6.8  ms"1  in  2000.  In  general, 
both  COARE  and  BVW  give  nearly  identical  results,  as  also 
evident  in  Figure  3b.  Annual  mean  difference  between  wind 
speeds  at  10  and  those  at  4  m  is  0.5  m  s'1  when  using  either 
COARE  or  BVW  in  adjusting  winds. 

5.  Differences  Between  Neutral  and  Stability- 

Dependent  Winds 

5.1.  Analyses  at  Buoy  Locations 

[26]  After  the  winds  are  adjusted  to  10  m,  as  explained  in 
section  4,  we  now  form  a  difference  between  stability- 
dependent  and  equivalent  neutral  winds.  This  is  first  done 
at  the  previously  analyzed  buoy  location  (0°N,  155°W)  in 
2000  and  2001  and  then  repeated  at  all  available  TAO, 
NDBC,  and  PIRATA  buoys  during  all  available  time 
periods  from  2000  through  2005. 

[27]  Equivalent  neutral  winds  are  calculated  using 
BVWN,  LKB,  and  LOG,  as  listed  in  Table  1.  Stability- 
dependent  winds  are  typically  weaker  than  their  counterpart 
neutral  winds,  which  is  consistent  with  the  typically  unsta¬ 
ble  stratification  at  this  buoy  (Figure  4).  Neutral  winds  from 
BVWN  and  LKB  are  almost  identical  in  these  examples, 
while  those  from  LOG  are  different  and  have  much  less 
variability  on  daily  timescales  in  both  years.  A  standard 
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Figure  3.  (a)  Time  series  of  daily  averaged  wind  speed  at 
4  m  above  the  sea  surface  as  obtained  from  a  TAO  buoy, 
(0°N,  155°W),  located  in  the  central  equatorial  Pacific  in 
2000.  Also  given  are  wind  speed  values  adjusted  to  10  m 
using  COARE  (version  3.0)  and  BVW  algorithms, 
(b)  Differences  between  wind  speeds  adjusted  to  10  m 
and  those  at  4  m.  Symbols  used  for  differences  (circles  and 
pluses)  are  generally  overlain  on  each  other  for  a  given  day, 
clearly  confirming  the  fact  that  both  algorithms  give  nearly 
identical  results  at  this  particular  TAO  buoy  location.  Note 
that  the  x  axis  is  labeled  starting  from  the  beginning  of  each 
month. 


LOG  correction  provides  a  neutral  wind  that  is  fairly 
consistently  <0.05  m  s_l  different  from  the  stability-depen- 
dent  wind. 

[28]  Hereinafter,  BVW  is  used  as  a  reference  algorithm 
since  it  has  an  option  of  calculating  either  stability-depen¬ 
dent  or  equivalent  neutral  winds.  Annual  mean  difference 
between  the  two  is  very  small  with  values  of  -0.09,  -0.07, 
and  -0.04  m  s_l  for  BVW  minus  BVWN,  BVW  minus 
LKB,  and  BVW  minus  LOG,  respectively,  during  2000. 
The  corresponding  differences  in  2001  are  again  small,  with 
values  of  —0.13,  —0.11,  and  —0.05  ms-1.  These  differ¬ 
ences  (annual  averages,  averaged  over  all  buoy  locations) 
are  small  for  most  climate-related  applications.  They  are 
substantially  smaller  than  the  accuracy  of  individual  wind 
speed  measurements,  0.3  m  s“ 1 ,  at  TAO  buoys  [Freitag  et  al. , 
2001];  however,  they  are  likely  to  be  statistically  significant. 
Regionally,  the  biases  can  be  larger  than  0.3  m  s  . 

[29]  Similar  to  the  TAO  buoy  location  examined  above, 
differences  in  stability-dependent  and  equivalent  neutral 
winds  are  also  analyzed  at  a  NDBC  location  (29°N, 
85°W).  This  buoy  is  west-northwest  of  Tampa,  Florida, 


United  States,  and  water  depth  is  «54  m.  Winds  at  this 
particular  buoy  location  are  measured  5  m  above  the  site 
level.  Hourly  measurements  of  near-surface  variables  are 
used  for  adjusting  winds  to  10  m  with  BVW,  and  for 
simplicity,  the  resulting  winds  are  then  averaged  over  a 
month  (Figure  5).  Also  included  are  monthly  averages  of 
equivalent  neutral  winds  computed  with  BVWN,  LKB,  and 
LOG  algorithms  to  show  differences  among  them.  Results 
are  given  from  2000  through  2005  so  that  one  could  see 
interannual  variations,  given  that  wind  speeds  and  stability 
can  change  depending  on  the  time  or  the  year.  At  this 
particular  location,  variations  in  wind  speeds  do  exist  from 
one  year  to  another,  but  the  common  feature  in  all  years  is 
that  winds  during  summer  are  typically  weaker  than  those 
during  other  time  periods. 

[30]  Differences  between  stability-dependent  and  equiva¬ 
lent  neutral  winds  calculated  with  BVWN,  LKB,  and  LOG 
are  given  in  Figure  6  at  (29°N,  85°W)  during  each  year, 
separately.  The  most  striking  feature  seen  in  Figure  6  is  that 
the  LOG  algorithm  is  the  outlier,  as  it  generally  yields 
neutral  winds  weaker  than  BVWN  and  LKB  values.  This 
difference  in  behavior  is  due  to  (1)  the  use  of  a  neutral  stress 
rather  than  the  stability-dependent  stress  and  (2)  the  fixed 
roughness  length.  The  differences  in  neutral  winds  obtained 
from  BVWN  and  LKB  are  relatively  small.  For  example, 
the  annual  mean  difference  for  BVW  minus  BVWN  (BVW 


Figure  4.  Differences  in  stability-dependent  and  equiva¬ 
lent  neutral  wind  speeds  at  a  TAO  buoy,  (0°N,  155°W). 
Comparisons  are  shown  for  (top)  2000  and  (bottom)  200 1 . 
In  all  cases,  stability-dependent  wind  speeds  are  adjusted  to 
10  m  (from  4  m)  by  including  the  effects  of  stability  using 
BVW,  while  equivalent  neutral  winds  were  adjusted  to  10  m 
using  BVWN,  LKB,  and  LOG  as  described  in  the  text.  Note 
that  unlike  the  others,  the  LOG  approach,  described  in 
section  2,  has  no  detailed  physics  in  it  and  uses  a  constant 
roughness  length. 
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Figure  5.  Monthly  mean  wind  speeds  at  10  m  above  the 
sea  surface  at  a  buoy  location,  (29°N,  85°W),  whose  NDBC 
station  ID  is  42036.  Stability  dependence  in  winds  is  taken 
into  account  using  near- sea  surface  atmospheric  variables 
in  the  BVW  algorithm  as  described  in  the  text.  Equivalent 
neutral  winds  are  computed  with  BVWN,  LKB,  and  LOG 
algorithms.  In  all  cases,  individual  hourly  data  are  used  to 
determine  wind  speeds  and  monthly  means  are  then 
constructed. 


minus  LKB)  is  generally  -0.2  m  s_l  (-0.1  m  s-1)  during 
most  of  the  years  from  2000  through  2005.  On  the  basis  of 
the  LOG  approach,  stability-dependent  winds  are  almost 
identical  to  equivalent  neutral  winds  over  the  annual  cycle. 
We  cannot  directly  conclude  which  algorithm  performs  the 
best,  but  on  the  basis  of  these  results  it  is  clear  that  BVWN 
and  LKB  give  comparable  results. 

[31]  Monthly  mean  differences  at  the  NDBC  buoy  loca¬ 
tion  of  (29°N,  85° W),  shown  in  Figure  6,  were  computed  on 
the  basis  of  hourly  wind  speed  measurements  as  mentioned 
above.  We  now  produce  a  box  plot  to  have  a  more 
qualitative  analysis  of  hourly  wind  speed  differences  used 
for  determining  monthly  mean  among  different  algorithms. 
As  an  example,  this  is  demonstrated  in  2001  (Figure  7).  The 
box  plot  simply  summarizes  hourly  wind  speed  differences 
at  the  NDBC  location.  Each  box  contains  the  middle  50% 
(i.e.,  median)  of  the  wind  speed  difference.  The  upper 
(lower)  hinge  of  the  box  indicates  the  95th  (5th)  percentile 
of  the  wind  speed  difference  and  is  limited  by  the  5th 
percentile  and  95th  percentile  differences,  and  the  whiskers 
show  the  largest  differences  within  the  month.  The  median 


Figure  6.  Differences  between  stability -dependent  and 
equivalent  monthly  mean  wind  speeds  at  10  m  above  the  sea 
surface  at  (29°N,  85°W).  Differences  are  based  on  values 
given  in  Figure  5. 


Figure  7.  Box  plots  of  differences  between  stability- 
dependent  and  equivalent  neutral  winds  at  (29°N,  85°W)  in 
2001.  All  these  are  calculated  using  hourly  stability- 
dependent  and  neutral  wind  speeds  for  each  month.  For 
example,  there  are  743,  670,  and  735  hourly  wind  speeds 
available  from  the  buoy  in  January,  February,  and  March, 
respectively.  These  winds  are  first  adjusted  to  10  m  using 
the  BVW  algorithm.  These  10  m  winds  are  then  converted 
to  equivalent  neutral  winds  using  BVWN,  LKB,  and  LOG. 
Differences  from  BVW  are  then  computed  to  determine  the 
statistical  summary  provided  in  the  box  plots. 
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Figure  8.  (a)  Time  series  of  monthly  mean  stability-dependent  wind  speed  values  and  equivalent  neutral 
wind  speeds,  both  adjusted  to  10  m,  from  various  places  over  the  global  ocean  in  2005.  Monthly  means 
were  formed  on  the  basis  of  daily  averages  at  each  buoy  location,  (b)  Differences  in  stability-dependent 
and  equivalent  neutral  wind  speeds  shown  in  Figure  8a. 


values  within  the  box  are  not  equidistant  from  the  hinges  for 
BVW  minus  BVWN  and  BVW  minus  LKB  during  most  of 
the  month.  Thus,  wind  speed  differences  are  skewed.  In  the 
case  of  BVW  minus  LOG,  median  values  are  very  close  to 
the  5th  percentiles.  Overall,  median  difference  values  are 
not  even  the  same  between  BVW  minus  BVWN  and  BVW 
minus  LKB,  but  they  are  very  close  to  each  other  (within 
0.1  m  s_I)  during  most  months.  This  indicates  that  LKB 
and  BVWN  generally  provide  identical  equivalent  neutral 
wind  speeds  even  though  they  have  different  physical 
parameterizations. 

[32]  To  give  a  further  idea  about  monthly  differences 
between  stability-dependent  and  neutral  winds,  more  exam¬ 
ples  are  provided  in  Figure  8  for  a  given  specific  year,  2005. 
These  locations  are  chosen  at  different  places  of  the  global 
ocean  so  that  we  can  examine  whether  or  not  differences 
can  be  quite  variable  from  one  location  to  another  since  the 
stability  conditions  may  vary  depending  on  the  region. 
These  six  particular  buoys  and  their  locations  are  as  follows: 
(0°N,  155°W),  a  TAO  buoy  located  in  the  central  equatorial 
Pacific;  (5°S,  165°E),  a  TAO  buoy  located  in  the  western 
equatorial  Pacific  warm  pool;  (12°N,  38°W),  a  PIRATA 
buoy  located  in  the  Atlantic;  (32°N,  75°W),  a  NDBC  buoy, 
whose  station  ID  is  41002,  located  in  the  east  of  Charleston, 
South  Carolina,  United  States;  (38°N,  123°W),  a  NDBC 
buoy,  whose  station  ID  is  46013,  located  north-northwest  of 
San  Francisco,  California,  United  States;  and  (58°N, 
150°W),  a  NDBC  buoy,  whose  station  ID  is  46080,  located 
in  the  Alaskan  coast. 


[33]  At  all  buoy  locations  we  first  adjust  hourly  wind 
speeds  to  10  m  using  BVW,  and  daily  means  are  formed. 
Equivalent  neutral  winds  are  computed  with  BVWN.  The 
conversion  to  10  m  buoy  winds  was  made  only  when  daily 
averaged  ocean  surface  temperature,  air  temperature,  and 
relative  humidity  were  available  from  the  buoy  measure¬ 
ments.  Otherwise  that  record  was  skipped.  Monthly  means 
were  formed  when  all  these  essential  variables  from  buoys 
were  available  at  least  20  d.  On  the  basis  of  some  tests,  a 
20-d  time  period  is  found  to  be  sufficient  to  represent 
monthly  mean. 

[34]  Figure  8a  shows  that  the  annual  cycles  of  equivalent 
neutral  winds  and  winds  including  the  effects  of  air-sea 
stability  are  typically  similar.  Equivalent  neutral  winds  are 
always  stronger  in  all  months  during  2005  except  at  (38°N, 
123°W),  where  the  stability  conditions  can  result  in  weaker 
equivalent  winds  in  May,  June,  July,  October,  and  Novem¬ 
ber  (Figure  8b).  Values  of  BVW  minus  BVWN  are  typically 
between  0  and  -0.3  ms-1. 

[35]  Finally,  we  would  like  to  examine  whether  or  not  the 
differences  between  two  wind  types  presented  at  six  buoy 
locations  are  typically  robust  at  other  buoy  locations  (see 
Figure  2).  Thus,  a  similar  procedure  for  calculating  stability- 
dependent  and  equivalent  winds  at  those  six  sample  buoy 
locations  (see  Figure  8a)  is  also  performed  at  all  TAO, 
NDBC,  and  PIRATA  buoys  from  2000  through  2005.  A 
scatter  diagram  reveals  that  there  are  no  systematic  differ¬ 
ences  between  the  two  (Figure  9a)  when  winds  are  adjusted 
to  10  m  using  BVW  and  equivalent  winds  are  calculated 
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Figure  9.  (a)  A  comparison  of  monthly  mean  equivalent 
neutral  versus  stability -dependent  wind  speed  based  on  all 
buoys  (NDBC,  TAO,  and  PIRATA)  during  2000-2005. 
Wind  speeds  at  10  m  are  calculated  on  the  basis  of  BVW 
and  BVWN  using  daily  wind  speed  values  at  each  buoy 
location.  Average  of  all  7619  monthly  mean  stability- 
dependent  (equivalent)  winds  is  6.3  m  s“r(6.5  ms"1).  (b)A 
comparison  of  BVW  versus  COARE  algorithm  based  on 
wind  speeds  adjusted  to  10  m  by  including  the  effects  of 
stratification.  As  in  Figure  9a,  monthly  mean  wind  speeds 
from  all  buoys  are  used.  Average  of  all  7619  monthly  mean 
stability-dependent  wind  speeds  adjusted  to  10  m  using 
BVW  (COARE)  from  all  buoys  is  6.3  m  s_1  (6.3  m  s-1). 


using  BVWN.  If  one  had  adjusted  winds  to  10  m  using 
COARE  instead  of  BVW,  the  same  results  would  have  been 
obtained  (Figure  9b),  indicating  that  the  use  of  either  one  of 
the  algorithms  in  adjusting  winds  does  not  change  results 
even  when  analyzing  results  at  all  buoy  locations.  This  also 
confirms  the  consistency  of  BVW  and  COARE  in  adjusting 
winds  to  10  m. 

[36]  Overall,  on  the  basis  of  76 1 9  monthly  mean  wind  speed 
values  at  all  buoy  locations  during  2000-2005,  neutral  winds 
are  stronger  by  0.2  m  s_1  than  stability-dependent  winds 
(Figure  9a).  As  to  wind  speed  adjustment  to  10  m  performed 


at  all  buoy  locations  during  2000  -  2005,  BVW  and  COARE 
yield  identical  results  with  no  (i.e.,  0  ms"1)  wind  speed  bias 
(Figure  9a).  This  also  explains  that  although  both  BVW  and 
COARE  essentially  involve  different  physical  approaches  (see 
section  2),  the  resulting  winds  adjusted  to  10  m  are  nearly 
identical.  Thus,  either  one  of  the  algorithms  can  be  safely  used 
for  the  adjustment  process.  Furthermore,  we  found  no  clear 
relationship  between  stability-dependent  winds  and  differen¬ 
ces  between  stability-dependent  and  equivalent  neutral  winds 
(Figure  1 0).  There  is  a  great  deal  of  scatter  from  a  wind  speed 
dependent  trend.  Thus,  at  a  given  location,  it  is  not  trivial  to 
estimate  the  difference  between  the  two  by  using  either 
stability-dependent  or  equivalent  neutral  winds. 

5.2.  Differences  Over  the  Global  Ocean 

[37]  We  now  examine  differences  between  stability- 
dependent  and  equivalent  neutral  winds  over  the  global 
ocean.  This  will  be  accomplished  using  wind  velocity 
measurements  from  the  SeaWinds  scatterometer  on  the 
QuikSCAT  satellite.  Hereinafter,  they  will  be  referred  to  as 
QSCAT  winds.  QSCAT  measurements  provide  an  excellent 
opportunity  to  examine  stability-dependent  and  equivalent 
neutral  wind  difference  because  (1)  they  are  calibrated  to  the 
neutral  stability  conditions,  thus  representing  equivalent 
neutral  winds,  and  (2)  they  can  be  converted  to  stability- 
dependent  winds  using  air-sea  stability.  Both  of  these  fea¬ 
tures  will  be  described  in  this  section  in  detail. 

[38]  Spacebome  radars  (scatterometers)  infer  surface 
winds  from  the  roughness  of  the  ocean  surface  [e.g.,  Liu , 
2002].  The  scatterometer  is  an  active  microwave  sensor  that 
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Figure  10.  Scatterplot  of  differences  between  stability- 
dependent  and  equivalent  neutral  10  m  wind  speeds  with 
respect  to  the  height-adjusted  stability-dependent  winds. 
This  plot  indicates  that  a  wind  speed -based  parameteriza¬ 
tion  cannot  be  used  to  approximate  the  differences  between 
stability-dependent  and  equivalent  neutral  winds.  Note  also 
that  most  low  wind  speeds  are  associated  with  stable 
stratification  and  that  the  influence  of  stratification 
decreases  as  wind  speed  increases. 
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samples  «90%  of  the  ice-free  ocean  in  1  d,  most  cells 
having  a  maximum  of  two  observations  per  25  x  25  km2 
grid  cell  each  day.  However,  there  is  a  latitudinal  depen¬ 
dence  on  sampling,  with  a  greater  average  at  higher  lat¬ 
itudes.  Measurements  of  radar  backscatter  from  a  given 
location  on  the  ocean  surface  are  obtained  from  multiple 
azimuth  angles  as  the  satellite  travels  along  its  orbit.  These 
observations  are  acquired  for  two  polarizations  at  different 
fixed  incidence  angles.  Wind  speed  and  direction  are  then 
inferred  from  measurement  of  microwave  backscattered 
power  from  these  multiple  look  angles.  QSCAT  was 
designed  to  measure  ocean  wind  vectors  with  the  accuracy 
requirement  of  root-mean-square  of  «2  m  s~ 1  for  the  wind 
range  of  3-20  m  s  .  These  requirements  are  essentially 
specified  for  a  cell  resolution  of  25  km  («0.25°). 

[39]  Remote  Sensing  Systems  provides  twice-daily  0.25°  x 
0.25°  fields  of  equivalent  neutral  winds,  available  online  at 
http://www.remss.com.  To  include  stability  effect  for  these 
winds,  we  first  interpolated  daily  averaged  sea-only  fields  (air 
temperature,  humidity,  etc.),  obtained  from  the  1.125°  x 
1.125°  resolution  European  Centre  for  Medium-Range 
Weather  Forecasts  40-year  reanalysis  (ERA-40),  described 
by  Kdllberg  et  al.  [2004],  to  the  0.25°  x  0.25°  QSCAT  grid. 
We  then  use  these  with  the  BVW  model  to  convert  rain-free 
equivalent  neutral  winds  to  stability-dependent  winds. 
ERA-40  is  used  when  available  (up  to  August  2002), 
and  the  1.0°  x  1.0°  Navy  Operational  Global  Atmospheric 
Prediction  System  is  used  afterward  (August  2002  through 
December  2005).  Note  that  sea-only  ERA-40  fields  were 
produced  by  only  using  ocean  points  (based  on  the  land-sea 
mask,  set  to  0.5  in  the  case  of  ERA-40 ’s  fractional  mask),  with 
creep-fill  interpolation  to  cover  land.  Kara  et  al  [2007b] 
provide  details  of  such  interpolation  to  construct  sea-only 
fields. 

[40]  The  process  of  forming  monthly  winds  is  identical 
for  equivalent  and  stability-dependent  winds.  In  our  pro¬ 
cessing,  only  twice-daily  rain-free  wind  measurements  from 
QSCAT  are  used  over  the  global  ocean.  First,  we  form 
monthly  averages  on  the  0.25°  x  0.25°  grid  using  a  cutoff 
of  20  rain-free  observations  per  month.  Then,  from  this  we 
produce  a  25-point  (1.25°  square)  observation-weighted 
average  at  each  0.25°  cell  using  a  cutoff  of  100  rain-free 
observations  per  month.  Finally,  we  fill  in  all  data  voids 
(land-  and  rain-contaminated  cells)  using  creep-fill  inter¬ 
polation  [Kara  et  al.,  2007b].  This  gives  us  a  data  set  on  a 
0.25°  grid  with  about  the  same  effective  resolution  as 
ERA-40  (1.125°  x  1.125°). 

[41]  As  an  example,  monthly  means  of  equivalent  neutral 
winds  processed  from  QSCAT  measurements  are  plotted  in 
Figure  1  la  over  the  global  ocean  during  January  2001 .  Note 
that  there  are  no  QSCAT  wind  observations  above  ice;  thus, 
regions  where  ice  is  present  (e.g.,  very  high  northern  and 
southern  latitudes)  are  masked  and  shown  in  gray.  The 
ice-free  regions  are  determined  from  an  ice-land  mask 
[Reynolds  et  al.,  2002].  The  mask  is  a  function  of  the  ice 
analysis  and  changes  by  month.  Equivalent  neutral  winds 
from  QSCAT  are  then  converted  to  winds  that  include 
stability  effects,  as  seen  in  Figure  11a  (middle). 

[42]  The  conversion  is  accomplished  using  6-hourly  near- 
surface  atmospheric  variables  (sea  surface  temperature,  air 
temperature,  and  relative  humidity)  obtained  from  ERA-40. 
For  the  atmospheric  variables,  we  first  formed  the  6-hourly 


fields  over  the  global  ocean  during  January  2001  and  then 
took  a  daily  running  mean  once  a  day.  The  same  process 
was  also  applied  to  other  months.  Spatial  variations  of  the 
resulting  stability-dependent  winds  are  generally  similar  to 
those  of  equivalent  neutral  winds.  This  is  also  true  for  the 
magnitude  of  winds.  In  fact,  the  difference  between  the  two 
is  very  small  over  most  of  the  global  ocean,  as  evident  in 
Figure  11a  (bottom). 

[43]  Similar  to  fields  produced  during  January  2001,  we 
also  form  monthly  mean  of  winds  during  July  2001 
(Figure  lib).  Again,  near-surface  atmospheric  variables 
from  ERA-40  are  used  for  determining  the  impact  of  air- 
sea  stratification  on  winds.  Seasonal  variation  in  wind 
speeds  is  clear.  For  example,  unlike  January  2001,  relatively 
strong  winds  are  not  seen  at  the  North  Atlantic  and  North 
Pacific  oceans  during  July  2001.  In  fact,  mean  monthly 
winds  are  generally  as  low  as  «7  m  s'1  in  these  regions. 
Spatial  variation  of  the  difference  field  between  the  two 
wind  types  is  also  quite  different  (Figure  lib),  although 
stability-dependent  winds  are  again  typically  weaker  than 
the  neutral  winds  except  at  very  high  northern  latitudes. 
Neutral  winds  being  stronger  by  >0.4  m  s_I  is  also  possible 
in  the  northernmost  Atlantic  and  Pacific  oceans.  As  in 
January  2001,  the  global  average  of  the  wind  speed  differ¬ 
ence  is  very  small,  with  a  value  of  -0.07  m  s-1  during  July 
2001,  but  spatial  variability  in  the  difference  values  must  be 
noted. 

[44]  An  investigation  of  the  differences  between  the  two 
wind  speeds  shown  for  the  months  of  January  and  July  in 
2001  is  extended  to  other  months  in  the  same  year  (Figure  1 2). 
The  common  feature  seen  in  all  plots  is  that  stability- 
dependent  winds  are  generally  weaker  than  equivalent  neutral 
winds  over  most  of  the  global  ocean.  However,  these  differ¬ 
ences  are  very  small  since  neutral  winds  are  typically  stronger 
by  0.1  m  s_I  than  stability-dependent  winds.  Spatial  patterns 
of  the  differences  are  generally  similar,  but  some  deviations 
arise  in  some  regions,  indicating  the  relatively  large  impact  of 
air-sea  stability. 

[45]  In  Figure  12,  differences  between  the  two  wind  types 
are  examined  only  during  2001 .  Are  differences  between  the 
two  similar  for  other  years?  To  answer  this  question,  we 
processed  QSCAT  winds  during  2000  and  2002.  Equivalent 
neutral  winds  were  converted  to  stability-dependent  winds 
using  near-surface  atmospheric  variables  from  ERA-40  as 
in  2001.  Figure  13  shows  differences  between  stability- 
dependent  and  neutral  wind  speeds.  The  resulting  differences 
are  again  small;  however,  there  are  regional  interannual 
differences  that  might  be  important. 

6.  Application  of  Stability-Dependent  Versus 
Neutral  Winds 

[46]  In  section  5.2,  we  found  that  stability-dependent  and 
equivalent  neutral  winds  are  typically  similar  to  each  other 
over  most  of  the  global  ocean.  Monthly  mean  differences 
are  generally  within  «0.2  m  s-1.  We  also  suggested  that 
stability  considerations  can  lead  to  important  differences  in 
wind  speed  on  short  timescales,  while  the  main  global 
results  are  presented  as  monthly  averages  only.  In  this 
section,  the  importance  of  air-sea  stability  effects  on  10  m 
winds  is  further  explored  with  an  application  involving  heat 
fluxes  over  the  global  ocean  on  shorter  timescales. 
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(a)  January  2001  (b)  July  2001 

Equivalent  neutral  wind  speed  Equivalent  neutral  wind  speed  (m  s-1) 


Stability-dependent  wind  speed  Stability-dependent  wind  speed 


Stability  dependent  minus  neutral 


Stability-dependent  minus  neutral 


Figure  11.  Comparisons  of  monthly  mean  equivalent  neutral  wind  speed  and  stability-dependent  wind 
speed  (both  are  at  10  m  above  the  sea  surface)  constructed  from  rain-free  QSCAT  measurements  over  the 
global  ocean:  (a)  January  2001  and  (b)  July  2001.  Note  that  monthly  means  were  constructed  from  daily 
fields. 


[47]  Our  main  goal  is  to  answer  how  air-sea  heat  fluxes 
would  change  if  one  used  equivalent  neutral  winds,  which 
are  directly  obtained  from  the  satellite,  instead  of  stability- 
dependent  winds.  As  an  example,  equivalent  neutral  winds 
from  QSCAT  winds  are  obtained  at  the  first  pass  of  the 
satellite  over  the  global  ocean  on  1  January  and  1  August 
2002.  They  are  then  converted  to  stability-dependent  winds 
on  the  basis  of  the  BVW  algorithm.  Unlike  the  monthly 
analyses,  differences  between  stability-dependent  and 
equivalent  neutral  winds  are  slightly  higher  (Figure  14). 
However,  the  largest  stability  effects  occur  in  ocean  regions 
that  are  characterized  by  either  western  boundary  currents  or 
meander  regions  of  the  Antarctic  Circumpolar  Current  (i.e., 
the  Southern  Ocean  regions)  as  noted  in  the  monthly  mean 
fields. 

[48]  Given  those  differences  in  stability-dependent  and 
equivalent  neutral  winds,  sensible  and  latent  heat  fluxes  are 
also  computed  on  the  basis  of  the  bulk  heat  parameterizations 
in  COARE  (version  3.0)  and  near-surface  atmospheric  var¬ 
iables  from  ERA-40  at  each  ocean  grid  point.  Heat  fluxes  are 
computed  with  two  inputs  separately:  (1)  stability-dependent 
winds  and  (2)  equivalent  neutral  winds.  Differences  in  fluxes 


are  formed  by  subtracting  equivalent  neutral  winds  from 
stability-dependent  winds.  The  resulting  sensible  and  latent 
heat  flux  differences  are  generally  negligible  (Figure  14). 
Sensible  (latent)  heat  flux  differences  are  within  1  W  m~2 
(2  W  m-2)  over  most  of  the  global  ocean. 

[49]  Including  the  impact  of  air-sea  stability  causes 
relatively  large  effects  on  heat  fluxes  only  near  the  Gulf 
Stream  and  Kuroshio  current  systems  on  1  January  2002 
(Figure  14a).  These  are  the  regions  where  winds  and 
magnitude  of  latent  heat  fluxes  are  relatively  large  in 
comparison  to  other  regions  (not  shown),  resulting  in  large 
differences.  Stratification  effects  on  winds  are  even  smaller, 
as  evident  from  the  heat  flux  difference  fields  on  1  August 
2002  (Figure  14b).  No  significant  spatial  variations  are 
noted  during  either  time  period.  Overall,  in  isolated  regions 
associated  with  large  air-sea  temperature  and  specific 
humidity  differences  (and  typically  ocean  fronts),  we  find 
that  the  adjusted  surface  turbulent  fluxes  differ  by  over 
5  W  m-2,  which  is  considered  significant  as  a  short-  to 
medium-term  bias  [  Webster  and  Lukas ,  1992].  These  biases 
vary  seasonally,  but  within  such  regions  they  are  often 
significant  as  judged  by  this  standard.  Other  work  in 


10  of  14 


C04009 


KARA  ET  AL.:  STABILITY  EFFECT  ON  WIND  SPEED 


C04009 


(a)  2001 

Stability-dependent  minus  neutral 


(b)  2001 
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Figure  12.  Differences  in  monthly  mean  10  m  wind  speed  (a)  from  February  through  June  2001  and 
(b)  from  August  through  December  2001.  Differences  in  January  and  July  are  given  in  Figure  11.  Blue 
(red)  denotes  regions  where  stability-dependent  wind  speed  is  weaker  (stronger)  than  equivalent  neutral 


wind  speeds. 

progress  shows  that  there  are  large  monthly  biases  in  these 
regions. 

7.  Conclusions 

[50]  Important  physical  processes  (e.g.,  surface  fluxes, 
wave  evolution,  and  in  some  cases  upper  ocean  mixing)  are 


typically  parameterized  in  terms  of  wind  speed  or  stress 
inferred  from  wind  speed.  For  historical  reasons,  wind 
speeds  are  used  rather  than  equivalent  neutral  wind  speeds, 
which  are  typically  available  from  satellites.  A  quantitative 
analysis  of  differences  between  the  two  wind  types  is 
desired  to  determine  possible  errors  (particularly  biases) 
due  to  improper  use  of  equivalent  neutral  winds. 
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(a)  2000  (b)  2002 
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Figure  13.  Differences  in  monthly  mean  10  m  wind  speed  during  February,  April,  June,  and  August  of 
(a)  2000  and  (b)  2002. 


[5 1 ]  Stability-related  differences  between  winds  and 
equivalent  neutral  winds  are  found  to  be  small  on  monthly 
timescales.  However,  sampling  is  sufficient  that  the  small 
signal  is  easily  identified.  When  scatterometer  winds  are 
used  for  the  calculation  of  surface  turbulent  fluxes  or 
surface  wave  models  or  for  product  validation,  they  are 
large  enough  to  contribute  to  physically  noticeable  differ¬ 
ences.  The  distribution  of  the  differences  between  the 
stability-dependent  winds  and  equivalent  neutral  winds, 
combined  with  the  number  of  samples  that  go  into  forming 
an  average  value,  reveals  typical  magnitudes  of  biases  in 
using  equivalent  neutral  winds  as  a  proxy  for  stability- 
dependent  winds.  These  biases  are  larger  on  short  timescales 
(e.g.,  hourly)  for  which  there  can  be  greater  departures  from 
neutral  stability  and  limited  sampling  of  independent  envi¬ 
ronmental  condition.  Furthermore,  conversion  of  equivalent 
neutral  winds  to  stability-dependent  winds  should  be  per¬ 
formed  with  a  technique  that  considers  the  dependence  of 


surface  roughness  on  stress.  A  simple  log  law  adjustment 
using  a  fixed  roughness  length  typically  estimates  relatively 
small  corrections  with  the  incorrect  sign.  The  LKB  and 
BVW  flux  models  have  very  different  physical  consider¬ 
ations;  however,  both  include  a  feedback  between  stability 
and  roughness  length.  The  differences  in  adjustment  be¬ 
tween  these  two  models  are  negligible,  typically  within 
±0.1  m  s"\  which  is  substantially  smaller  than  most 
adjustments. 

[52]  Wide-ranging  applications  of  QSCAT  winds  would 
benefit  from  the  conversion  of  equivalent  neutral  winds  to 
stability-dependent  winds  which  are  used  for  computing 
wind  stresses.  For  example,  winds  are  generally  weaker  than 
the  neutral  winds  by  ±0.2  m  s“ 1  over  the  global  ocean,  but 
differences  between  the  two  can  even  be  up  to  ±0.5  ms-1  at 
high  latitudes  on  a  monthly  timescale  and  much  larger  on  a 
shorter  timescale  (e.g.,  daily).  For  comparisons  between 
QSCAT  winds  and  model  or  buoy  winds,  on  a  shorter  scale 
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Stability-dependent  minus  neutral 


4  -0.4  43  -0 2  41  0.0 


0.1  U  0.3  0.4  4 


Sensible  heat  flux  difference 


<  4  3  2  1  0 


2  3  4  4 


4  4  3  2  1 


0 


2  3  4  4 


Latent  heat  flux  difference 


4  -8  *6  4  *2  0 


Latent  heat  flux  difference 


4  -8  *6  -4  *2  0  2  4  6  8  4 


Figure  14.  (a)  Differences  in  10  m  wind  speed  and  sensible  and  latent  heat  fluxes  during  the  first 
satellite  pass  on  1  January  2002.  Differences  in  heat  fluxes  are  computed  using  stability-dependent  and 
equivalent  neutral  winds  in  the  CO  ARE  (version  3.0)  algorithm  as  explained  in  the  text,  (b)  The  same  as 
Figure  14a  but  during  the  first  satellite  pass  on  1  August  2002. 


(e.g.,  hourly  or  daily),  there  are  likely  to  be  sufficient 
observations  that  ignoring  this  adjustment  will  result  in  a 
statistically  and  physically  significant  bias.  Our  previous 
work  has  also  revealed  that  biases  due  to  ocean  current  and 
wave  motions  can  be  regionally  important  [Kara  et  ai, 
2007a]. 

[53]  For  potential  users  of  interest,  hourly  and  daily  wind 
speeds  (both  stability-dependent  and  neutral  ones)  and  other 
near-surface  atmospheric  variables  adjusted  to  10  m  using 
historical  data  (1980  through  2005)  from  all  TAO,  NDBC, 
and  PIRATA  buoys  are  available  (http://www7320.nrlssc. 
navy.mil/nasec).  The  adjustment  is  made  using  all  air-sea  flux 
algorithms  (COARE,  LKB,  BVW,  BVWN,  and  LOG)  dis¬ 


cussed  in  this  paper.  Rain-free  QSCAT  winds  are  also 
available  over  the  global  ocean. 
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